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1. llie area of research of this grant vas the study of a novel 
nethod of extracting sea height information - the sanple functions, 
generally - from SAR conplex data, a method that vas suggested by a fun- 
damental SAR ocean imaging model for gravity vavesll], that showed that 
information about the long vave is present in the SAR conplex data, espe- 
cially its phase. The initial study [2], [3], mainly supported by earlier 
NASA and (X® grants and partially supported at the publication stage by 
the present grant, enployed an ad hoc, 2 dbeit quite reasonable, phase 
demodulation algorithm followed by linear regression amd filtering: only 
the latter two steps incorporated a priori information that might be 
available. A relatively siirplified simulation indicated that, among other 
possible limitations, the finite bandwidth of the SAR system inposed the 
apparently most serious limitation. A preliminary application to 
SEASAT-SAR cotplex imagery was encouraging. This work has been reported 
to the oonmunity, formally and informally [2], [3], [4], [6]. It was evident 
that a priori information about the long wave, at a minimum, if incor- 
porated into a "more sophisticated" piuise demodulation structure - i.e. , 
at an earlier point in the algorithm, could conceivably mitigate this 
bandwidth limitation. 

2. The work on this grant began by directly addressing the following 
central problem: given the SAR cotplex data of the sea, modeled as 
described above, and received along with thermal noise, v/nat is the opti- 
mal (minimum msan-square error) estimator (a conditional expectation, 
then) of the long wave structure and vhat is its performance? While such 
an estimator is generally very difficult to find and, when found, imple- 
ment, nevertheless ve have done so for this problem! The significance of 


sudi an answer is sinply that it gives the best possible performance wha- 
tever deleterious effects may be modelled - finite bandwidth, thermal 
noise, random phase, etc. 

(a) In particular, it has been found that this c^timal estimator is 
able to overcome, to a considerable degree, the bandwidth limitation 
encountered by the earlier-posed, ad hoc, sub-c 5 )timal estimator. The 
optimal estimator is able to do this because it incorporates in its struc- 
ture the a priori knowledge available). This estimation structure can be 
quite flexibly adaptive (at the cost of increased ccxrputation), is 
realized as an efficient, recursive calculation - e.g., as inferred from 
the SAR image, and, generally speaking, seems practical. 

While we will defer detailed discussion tu a forthcoming article [8], 
we call attention to the simulation shown in Fig. 1. The complex image 
was (randomly) generated using our model of the SAR ocmplex image of the 
sea and was processed the cptimal estimation structure to produce an 
estimate of the height of a sinusoidal long wave - of known jAiase and 
wavenumber, as could be separately estimatd from the SAR image, e.g. the 
long wave shape is shown in Figure 1 and the resulting mean-square error, 
and its Sctmple average, for 50 simulation runs is shown in Figure 2. The 
wavenumber spectrum generated by the SAR ocean-sensing mechanism is 
several times greater than the SAR bandwidth but the performance is very 
good! ! 

(Also shown for reference is the Cramer-Rao lower bound when the height is 
an unJaxuwn parameter: in the simulation it was assumed to be a Rayleigh 
random variable with an a priori known variance. ) 

(b) The above-discussed (^imal model has, in the main, been ela- 
borated so far atssuming a stationary scene. Actually, it may well turn 
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exit that the assurrption vas not especially limiting in view of analysis 
concerning the image model nature. Under recisonable conditions (e.g., an 
I/-band SAR with not too fine a resol\ition and typical sea parameters) 

(i) the small wave structure part that can influence the SAR image is more 
concerted than dispersive in its action, and (ii) an appropriate focus 
adjustment cem "render the long wave stationary". Then, with’/i some 
details, the SAR imaging model reduces to that used in the above-discussed 
height estimation study. 

(The observation (ii) is a well-Jowwn controversy in this comnonity 
euid ’..e have had a fairly general proof that the focus adjustment is deter- 
mined by the long wave's phase velocity for scxne time. The ejqjerimental 
"test" of this focus peurauneter dependence, proposed at an APL workshop 
last October was, as reported at a fftRSEN worJeshe^ at JPL in January, in 
each instance supportive of our SAR model: that is, that the dependence is 
on the long wave's phase velocity. ) 

(c) In considering the practical inplementation of this optimal 
estimator, it was noted that a significant saving in oonpitation can be 
achieved by generating the short wave ensanble by a so-called "diaotic 
dynamical" process - rather than the "conventional" (e.g., Markov) random 
process models. Such am observation is of much wider consequence and em 
initial publication has been prepared [7]. Besides its potential com- 
putational advantage [8], such "chaotic'- models of the sea surface are 
known to arise naturally as solutions to the nonlinear hydrodynamical 
equations and, hence, may be precisely the kind of models needed. 

(d) The nature of the short wave ensemble plays a critical role; 
its presence is necessary for receipt of backscattered energy at typical 
intermediate incidence angles; at small incidence euigles a quasi-specular 
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bbi?kscatter from the long vave can be more significant. The structure of 
optimal height estimator, and its performance depends upon the sta- 
tistical nature - e.g., "odierence length" - of the short vave ensemble 
and on the backscatter mode eis determined by the incidence zuigle. A simu- 
lation of the estimation algorithm - and the SAR image - in a study vgould 
be very Informative, as vould be eventual ooirparison with the data forth- 
coming from the SIR-B e^ 5 >eriment, offering data at various incidence 
angles. 

(e) The algorithm was simulated through numerical simulation of the 
optimal estination algorithm to estsablish, as conpletely as possible, its 
accuracy, flexibility, and practicality; The attenpt to apply the 
algorithm to SEASAT-SAR data, as supplied by NASA JFL and/or ERIM, was 
only partially acconplished: despite repeated requests to both, only a 
limited amount of data was obtained from ERIM: see [8], included in this 
report. 

3. Hie following publications, presentations, and discussions - 
acconplished and prospective - have been done during the Grant period. 

(It is likely that the topic of [8] will produce several publications - 
21 S, e.g., [7].) 

I. Journal, book, and proceedings publications supported by the Grant: 

(1) "The SAR image of short gravity waves on a long gravity wave", 
in Proceedings of a Syitposium on Vfeve Dynamics and Radio 
Probing of the Ocean Surface, O. M. Phillips £uid K. Hasselman, 
Eds., Plenum Press (in press). (Partial support of revision 
and manuscript pr^>aration; also supported by earlier NASA 
grant . ) 
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(2) "A sea surface height estimator using synthetic aperture radar 
complex imagery", IEEE Trans. Ooean Engineering, ^ril 1983. 
(Partial support for revision and manuscript preparation; also 
supported by earlier NASA grant. ) 

(3) "A sea surface height estimator using SAR conplex imagery". 
Proceedings of Oceans *82 Conference, Vbshington, DC September 
1982. (Partial support for travel and manuscript preparation; 
also supported by earlier NASA grant.) 

II. Hbrkshop participation related to grant: 

(4) "A novel SAR spectral estimation algorithm", presented at 
SEASAT-SAR Wbr)cshop CXI Ocean Vbve Spectra, Johns H^)kins Applied 
Physics Laboratory, Oc±ober, 1982. 

(5) ONR Workshop cxi SAR Oc^n Imaging Applications, Jc^ns Hopkins 
Applied I^ysics Laboratory, Decretiter, 1982. 

(6) NASA LARSEN VforkshcTp on SAR Ooean Iiteging Theory and Ejq^eriment, 
Cal. Tech. JPL, Jan., 1983. 

(7) "Optimal estimation vd.th chaotic dynamics", Prcxreedings of the 
1983 Conference on Information Sciences and Systems, Jcrfins 
Hcpkins University, March, 1983. 

III. Publications in prcxass: 

(8) "A Fundeunental Model and Efficient Inference for SAR Ocaam 
Imagery", being revised for publication in IEEE Jo. Ocaanic 
Engrg. 

Enclcaed are a reprint of [7] and a preprint of [8]. 
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SUM1ARY 

Recently deepening understanding of nonlinear dynamical systems has 
revealed an interesting aspect of their possibly oonplicated behavior - 
their so-called 'chaotic', seemingly random, evolution. The observation 
of such an evolution leads naturally to the problem of estimating aspects 
of such systems, a problem area of increasing inportance as the use of 
such models spreads into many sciences. 

While the evolution of diaotic systems can be very complex, it is the 
result of Iterating a nonrandom mapping of a (possibly random) initial 
state: optimal estimation structures can thereby simplify, relative to 
conventional reuxiom evolutions, possibly to practicality in specific 
instances. 

This laist remark is of additional significance since, in practice, it 
is often true that only qpiite limited knowledge is available for model 
construction - e.g., a correlation function: the 'orbits' of a chaotic 
dynamical model may well have a suitable sanple-path correlation. 

Typically, the evolution of such chaotic dynamical systems is indexed 
by a 'chaos parameter': both the gross and detailed nature of the evolu- 
tion can depend upon its precise value. Thus a oonputationally simple 
dynamical model is capable, under parametric control, of generating a 
diversity of 'sample functions'. Of course, the obverse is that, first. 
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care oust be taten in applying sudi nodela and, second, the performance of 
estimstion sch^nes nay also sensitively d^nd upon the value of the chaos 
parainetcir. 

The fomulation of a recursive estimation problem, in discrete tine, 
using Markov and chaotic descriptions is contrasted. Iben a sinple para- 
meter er<tiination problem is noted in whid) the more general role of cer- 
tain orbit functionals appears. Next a more conplicated nonlinear 
filtering problem - arising, e.g., in an ocean ronote sensing system - is 
discussed. Finally, the 'identification* of the chaos parameter itself is 
studied, assuming an a priori distribution. G^^nerally conditional espec- 
tation estimators are considered vdth the assistance of numerical 
oonputaticMi. 
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Abatract 

Nonlloaar aapploga that can axhlbit 
“chaotic.'* aeemingly random, evolution have 
appeal aa modela of dynamical ay sterna* Their 
deu^mlnlatic evolution, via-a'^vis Markov evolu- 
tiona, can result in much simpler optimal 
detection and estimation algorithms. The varia- 
‘ tion of a “chaos” parameter (m) can result in 
I diverse evolutions, suggesting a simple but rich 
source of model variations. For the specific 
mapping examined here, this latter possibility 
is problematic due to the extreme sensitivity 
•on u of the evolution in the "chaotic regime." 

1. Introduction 

Recent deepening understanding of nonlinear 
systems has revealed an Interesting aspect of 
their possibly complicated behavior - their so- 
called "chaotic.” seerlngly random, evolution. 

The observation of such an evolution leads 
naturally to the problem of making inferences 
regarding aspects of such evolutions, a problem 
increasing in importance with the spreading 
study of such models in many scientific areas. 

Uhile the evolution of such a system can be 
very complicated indeed in detail, in one respect 
it is simple, evolving as an iterated, nonrandom 
mapping on a given initial state. Optlsul algor** 
ithms for inference - detection and estimation - 
can be computationally quite feasible. 

This observation and the remark that, in 
practice it is often true chat only quite limited 
knowledge - e.g.. a covariance function - is 
availsDlc for model construction, leads to the 
idea that a nonlinear, "chaotic” model might 
provide evolutions with suitable properties. 

Such chaotic dynamical models may be indexed 
by a "chaos" parameter: both the gross and 
detailed nature of the evolution may depend 
dramatically upon its precise numerical value. 
Thus a relatively simple dynamical model can be 
capable of generating a diversity of evolutions, 
under a simple parametric control. Paradoxical^, 
j careful consideration would be required in a 
ipractical application. Also, the performance of 
■optimal Inference algorithms may depend sensitiv- 
ely on this parameter. Both matters will be 
addressed here. 

2. Harkov and Chaotic Inference Models 

(a) Consider a discrete time, nonlinear 
system evolving in accordance with the dynamical 
model 

I x(k+l)-f[x(k), «(k)l, k-0.1,;.., (1) 
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where {w(k)} is a ssquenceof mutually indepen- 
dent random variables ("r.v.s”). The evolution 
(1) is observed ss 

yC.)-£lx(k).v(k)]. (2) 

%fherc (v(k)) is s sequence of mutually independ- 
ent r.v.s, the random processes iw(k)} and 
(v(k)} being independent, (f and g are real 
functions of two real arguments.) 

Taking up, e.g., the filtering problem - 
optimal estimation of x(k) given Che observation 

y(k) ; (y(0) y(k)). the knowledge sufficient 

for estimating x(k) in accordance with common 
optimality crlterla-e.g. , minimizing the estim- 
ate's mean-square error ("MSE") - is given by 
j^he a posteirorl distribution of x(k) given 
y(k), denoted p(x(k)| y(k)]. (We assume density 
functions ("d.f."s) exist in every case.) It is 
well-known that the evolution of this a poster- 
iori d.f. is governed by the recursion 

r /f^(k+l ,k)dx(k) . (3a) 

o[x(k+l)ly(k+l)I- ;p(|,+i^it)dx(k)dx(k+l) 

where 

F(k+l.k)=ply(k+l)fx(k+l)) p(x(k-fl) |x;k)l 

p(x(k)|y(k)l. (3b) 

It is also well-known that this recursion is. 
generally, very difficult to specify in detail 
and implement as a feasible computation. The 
jinimum MSE estimator ("M-MSE-E") of x(k), given 
y(k), is then 

i(k(k) . /x(k)p[x(k)i;(k)ldx(k) (3c) 

where P[x(k)|y(k>] is computed recursively by 
(3a). 

(b) Alternatively, consider the evolution of 
a possibly "chaotic" nonlinear system defined by 

x(k+l)-f ^(x(k) ] ,k.0,l. . . . ,x(0)-Xq. (4) 

where is an initial state and p is a ("chaos") 
parameter; Xq and/or p may be random variable*^. 
("r.v*"s). Conditioned upon any such randomness, 
if present. (4) describes a non^-andom evolution 
that, when replacing (1) as a dynamical mode], 
results in a great simplification of (3), namtjly 

p(x(k+l)|y(k+l)] - (5a) 

//p(y(k+l)|x(k+l)l6(x(k+l)-f‘^*^(xQ))C(k)dXQdi) 
;/G(k + l):dxQdp 
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(5b) 


{ C(b)iply(k) l*(k)-f|j(xQ) ]C(k-l) .C(0 )*p(xq)p(m) . 

' Here x(k/-* 

fO(x )=x ^ ® 

i *0‘ 

I Thus. s.g.. the M-MSE-E of x(k)» givsn y(k). 

. <5c) 

i(k 1 k)-//f “ (»Q)C(k) d»pdu///C (k)d»(,di, . 


Here x(k)«f (x^) Is the k-th iterate of f . 
" M 0 M 


A variety of related inference probleas can 
be posed, all aiaplifying greatly as a consequ^ 
cnce of the siapllficstion of the joint d.f.: 


; p[y(k),x{kkul»» Ply(i)lx(l)l . 
i ^ 1-0 

1 » «{x(j>-fj(k*)}' p(Xq)p(w). 

I J »1 . ^ “ 


(5d) 


(c) The model <4) can result in far less 

I calculation than the aodcl (1). Suppose x(k |k) 

I at k-k* is required: (3c) requires Ak* Integra- 
! cion; in contrast. (5c) requires only 4. If 
x(k) in model (1) is s vector of dimension N» 
the integrations required are. reap.. 4k*N and 4. 
This latter comparison is of interest if. in 
such a case. (4) can. n sc7»e sense, adequately 
replace (1) as a model. 

(d) For the specific numerical calculations 
here we choose the mapping 


f^(x)sl-px^ , xe[-l.ll,MEl0,2) ; (6) 


it is representative of a well-studied class of 
mappings that are (1) continuous, (ii) of one 
maximum (at x«0)» (iii) monotone decreasing with 
: increasing |x|» and (iv) of a certain convexity 
in the derivative 3f|3x[ll. The general nature 
of the "orbits" {x(k)-f“(xQ) , k>0) is suggested 
by Fig. 1 where {x(k). k*500 to 600) is plotted 
versus u. Viewing it. scxse of the following 
facts are agreeable. 


I 
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Fig. 1-The iterates (500 to 600) of (6). 
For sufj^iclently small values of y. Che 


orbits, for almost all xq. approach a "final" 
value x("). a fixed point of the mapping f^. As 
y is increased sufficiently, e "bifurcetion" 
occurs, the orbit then, for almost sll xq, 
approaching a stable, periodic visit to two 
values, fixed points of the mapping fp. As y 
increases, further bifurcations occur on an 
increasingly finer scale, until y ■ 1.40155.^. 
ia reachtd. where a quite different behavior is 
encountered. 

There an aperiodic motion on a set of 
(Labaaqut) measure 0 occurs; the orbits of almoft 
all xq arc attracted to this act. and the orbits 
initially "close" remain close (for almost all 
Xq) - an ergodic. but not mixing, evolution. 

When y«2, the orbits arc also almost all 
speriodic, but range over the entire Interval 
[-1.1): an invariant measure, absolutely contin- 
uous with respect to L.-.hcsque measure, exists 
and its form is known [2]. Orbits initially 
close almost surely do not remain close-an 
ergodic and mixing evolution. (Incidentally, the 
points XQ»0.il/2.±l are easily directly seen to 
be exceptional polnts-as a numerical analysis 
may well Inadvertently discover!) 

It is not known if other types of evolution 
occur: in fact* the dynamics of this model are 
not fully understood, though s great deal is 
known {!)• The following facts are helpful in 
orienting and evaluating numerical analysis 
(i) f has either one or no stable, periodic 
orbir . (ii) If f has a stable, periodic 
oribit. then the orbits of almost all xq are 
attracted to It and, specifically, xq-0 is so 
attracted. 

3. A simple estimation problem 

A typical question that arises concerning 
the nonlinear model (4) occurs in the following 
simple estimation problem. Suppose that xq md 
y are known s priori and that (2) Is specialized 
to 


y(k)-a x(k)*h^(k), k-0,1,,.., (7) 

where a Is an unknown parameter to be estimated, 
having data y(k). Suppose further that the 
(v(k)) are identically distributed normal r.v.s. 
with mean 0 and variance . 

Then it is easy to show that the ouixlmum 
likelihood estimator, 

S(k)-( I y(l)t^(X(,)l/i r lfj(x )]^), (8a) 

1-0 i-0 “ 

is efficient- that is, it is unbiased and of 
error variance equal to the Cremer-Rao lower 
bound, namely 

E{(l(k)-.)^)- (8b) 

Therefore, of special interest is the average 
along the orbit, ^ 

V(k;x ,w)i ^ 

1-0 
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Cirtaln V»Mtions arit«: ••g*$ fot 

larf« k« do«a V(k;xQ»y) dapcnd upon and/or y? 
' For cha nodal 16K V(k;xQ,y) la aaaily nun- 
arically calculacad. Tha raault can ba axpactad 
to ba influancad by two phanonana: first, a 
nuBbar of Itarations nay ba raqulrad bafora tha 
orbit arrivas on tha "attractor" (assuning it 
axiata) and, aacond, a vary large nunbar of itar- 
jations nay ba required to achieve "atabla" 
'average value over the attractor* Ilia first 
.affect is noticeable but hare nlainisad by sunn- 
iog fron tha 100th iteration on. A sun over tha 
jaubaaquant 100 itarations resulted in an orbit 


javaraga V apparently independent of xq for trial 
Ivaluaa of y ■*4,.8, and 1.2: sea Fig* 2 where 
|xq>*(0,0.05, . . . ,1). Further calcula.ion over 
{1,000 (vice lOO) itarationa resulted in indepen- 


dance of xq of V whan y«l*6 and, whan y«*2, furdiar 
calculation over 5,00u itarations raaultad in a 
stable value of V, independent of xq. 


ware (0.3a,0.332....), (1,1. 01*..) (l.<»6, 1*45.*.). 
(1.76,1.763..*), all for v-0*l. These Simula- 
tlons enployad a net of 100 points in y: as k 
increased, tha aupport uf the a posteriori d.f. 
p[y|y(k)] became . within the numerical range of 
tha computer ( 10i43S), confined to a small 

portion of these points; *V , a more refined, or 
adaptively refined, y-net would be of interest. 


Fig* 3- Tha H-MSE-E y of (9) for several y. 


i 

t 

i 

?ig*2-The orbit average V of (8c) versus x^, with 
y as paranetar 


The general conclusion at this point is chat 
the parameter y is efficiently and accurately 
estimated by the M-MSE-E. At issue, however, 
is a aomewhat more subtle matter, as will be seen. 

Cramer-Rao bound . - The Cramer-Rao lower 
bound on the naan-square error that any estimator 
may have, is of interest-chough it is not nec- 
essarily attainable by any estimator. This 
lower bound involve the generally well-known 
form-cake Xq and y to be unknown parameters - 
2 k 

-E( ^ I 

3 / 1-0 ^ ^ 


4* Chacs Parameter Estimation 

' Taking up the "chaotic" dynamical model (4) 
and the nbservatlon node! (6), (7) - with a*l, the 
M-MSE-E y(k) of y, given data ^(k), is, by (5), 


Here, given (6), the derivative with respect 
to y can be defined recursively as 

(10b) 


^ f^ 
3y 


«y 


-1 3 f°-0. 


y(k)« /yp(y(k) |y]p(y)dy//p(f(k) |y]p(y)dy (9a) 


where 


(9b) 


p(y(k)liil- 


✓21,, 


I (y(i)-fJ(K )1^). 
2v i-0 


The euaa are efficiently, recursively calculated. 

The dynamical, observation; and optimal 
estimator equations were numerically simulated 
and some results are shown in Fi^. 3, for y c 
(0.33.1,1.46,1*76), all with x -0. Referring to 
Fig* 1, the rate of convergence for y e (0*33,1, 
|l*76) intuitively corresponds to the "dlstlnct- 
iiveness" of these y values; the quite rapid 
jconvaigence for y*l*46 in the chaotic region is 
informative. . The (true value, estliaate) -pairs 


The behavior of 3f^/3y itself is of some intecust: 
like f^, it displays a rich behavior* While 
bounded-lt is a polynomial in y of degree less 
chan n, it can be relative^ large at points 
(in y) of bifurcation of f , and it is very large 
in the "chaotic region." In Fig. 4, the 90th 
through 100th iterates of this derivative are 
plotted over a restricted y-set in the "chaotic" 
regime so that detail may be seen. It can be 
many orders of magnitude greater still at larger 
y. At least some of the nature of 3f”/3y is 
inferrable from the nature of as evidenced 
in Fig. 1* 
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corrclAtlon function - would then be known, the 
conclusion sppesrs to be unjustified! To cler- 
ify this, sharper Mssurea of “orbit dlBtingulah> 
ability*' are now considered. 

5. Orbit distin&uishability 

(a) The results of the chaos paraacter 
estiaation problem Just discussed allowed the 
possibility that the orbits of (6), as observed 
via (7) (with a*l), arc extrenely distinguishable. 
To aharpan this question, suppose one of two 
orbits, corresponding to and U2* accord* 
ingly observed: then, by any of the usual decis- 
ion criteria, the decision perforaance is deter- 
ained by the “dietance" 


Fig. 4 - The derivative df^/dy of (lOb) , for n • 
to 100, versus m* 

The orbit average appearing in (10a), 

is shown in Fig. S for k«20, x^-O.OS, and 40 
points in M. As the nuaber of^y points are 
: increased, greater detail in <k>'a dependance un 
;y is seen in the chaotic region, generally to the 
llait of the computer graphics and beyond. The 
nature of the y-dependence otherwise is intuit- 
ively agreeable with the dependence exhibited in 
Figs. 1 and 4. 



2 

noraalised by v . 

In Fig. 6 is shown this distance, depending 
on a sequence yjj = y.-(0.1)^ with y»1.6 and for an 
orbit segment from^k-lOO to 1,100. The remark- 
able result is chat, to the accuracy of the 
computer's representation of numbers (12 signif- 
icant digits), the orbits remain distinguishable! 


Fig.5-The orbit average (lOc) appearing in the 
Creamer-Rao bound. 

The very large values of <k>, for this rel- 
atively small k, for most of the y in the chaotic 
region, seem to there allow extremely accurate 
estlsation of y. (A Monte Carlo simulation of the 
M-HSE-E (9) could establish, numerically, whether 
or not such estimators exist.) But such a con- 
clusion would depend upon the utilization of the 
estimate: if it Is to be able to “identify", or 
fix, the chaotic dynamical system so that the 
lorbit, and itp attendant proper Li£& - e.g. , orbit 


Fig. 6- The distance (11) between orbits of f 
snd fpj. y 2 -y-( 0 . 1 )I. ^ 

(b) Another TOasure of distinguishability 
of the orbits is the “orbit correlation coeffic- 
ient" 



Referring to Fig. 7, it is seen chat this 
same extreme distinguishability persists to the 
limit of machine accuracy. (The values of unity 
for I near 12 occur when the machine can no 
longer distinguish p.and y 2 *) At the same time 
the correlation wlth^orbits considerably removed 
in y persists. (As a check on the hoped-for 
Independence, of the calculation on xq, the curve 



« di»c«nc« btcw«(in o.c.f*M* ratlmr thani uu 
•arlitr* • di«tanc« bttwcen orbica themaelves. 
Writing this fora out» 


7 t I X rk/ \ 

.Ll'^ „"h '■»’ 




In Fig. 8 this distance is shown for s 
sequsnea (I/IO)^ approaching 

the “chaotic region. “ for *^Bin"^®0» 

L*2S« Again, to the computer's limit of 

12 significant digits, the orbit spectral densi- 
ties ar^- distinguishable t 


I 


i 


Fig. 8- The distance (13) between orbit spectral 
densities of f,, and f for w_» 


(d) This extreme sensitivity of the orbit 
of (6)i and its orbit correlation function, to 
the precise value of the parameter p in the 
chaotic region renders probles&atic the feasibil- 
ity of using such a model in a practical case. 
Indeed, it raises a serious question whether 
the oft-cited and much-studied nonlinear mapping 
(6) is useful as a physical model at all. 
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